Reheating Effects in the Matter Power Spectrum and Implications for Substructure 
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The thermal and expansion history of the Universe before big bang nucleosynthesis is unknown. 
We investigate the evolution of cosmological perturbations through the transition from an early 
matter era to radiation domination. We treat reheating as the perturbative decay of an oscillating 
scalar field into relativistic plasma and cold dark matter. After reheating, we find that subhorizon 
perturbations in the decay-produced dark matter density are significantly enhanced, while subhori- 
zon radiation perturbations are instead suppressed. If dark matter originates in the radiation bath 
after reheating, this suppression may be the primary cutoff in the matter power spectrum. Con- 
versely, for dark matter produced nonthermally from scalar decay, enhanced perturbations can drive 
structure formation during the cosmic dark ages and dramatically increase the abundance of com- 
pact substructures. For low reheat temperatures, we find that as much as 50% of all dark matter is 
in microhalos with M > 0.1M® at z — 100, compared to a fraction of ~ 10 -10 in the standard case. 
In this scenario, ultradense substructures may constitute a large fraction of dark matter in galaxies 
today. 



I. INTRODUCTION 

The standard cosmological model accounts for a re- 
markable span in the Universe's history, from the epoch 
of big bang nucleosynthesis (BBN) to the era of struc- 
ture formation to our present dark-energy driven abyss. 
In the standard paradigm, inflation [lH3[ cures the hori- 
zon and flatness problems of the hot big bang and gen- 
erates the spectrum of curvature perturbations needed 
for the formation of large-scale structure. Inflation ends 
when the kinetic energy of the inflaton field is no longer 
small compared to its potential energy. At some point 
after the end of inflation, the inflaton decays, and the 
Universe becomes radiation dominated. 

The inflationary paradigm does not uniquely specify 
the reheating and thermal history of the Universe, how- 
ever. The only constraint on the Universe's temperature 
at the beginning of the radiation-dominated era comes 
from the thermal production of neutrinos; lowering this 
temperature lowers the neutrino abundance 0-0] > which 
has ramifications for BBN, the cosmic microwave back- 
ground (CMB), and large-scale structure. If the Universe 
is radiation dominated at a temperature of ~ 3 MeV, then 
thermal neutrino production is sufficient to produce the 
observed abundances of light elements [4|47| and the ob- 
served CMB and matter power spectra (8, 9], provided 
that the Universe is baryon-asymmetric at this temper- 
ature with a baryon-to-photon ratio of r\ ~ 6 x 10 -10 
0, E3- While challenging, it is possible to generate 
this baryon asymmetry if the Universe is not radiation 
dominated until it reaches a temperature of ~ 10 MeV 
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(e.g., [Ill H2|). 

At the end of inflation it is usually assumed that the 
inflaton can be described as a coherently oscillating field 
that eventually decays and transfers its energy to a ra- 
diative plasma [Ia4l9j . The inflaton is the first dom- 
inant energy component of the Universe and its decay 
could be the start of the standard hot big bang — or this 
could be a transient epoch. It is possible that a specta- 
tor field during inflation comes to dominate the energy 
density of the Universe after the inflaton decays; this 
scenario is a generic consequence of stabilized moduli in 
string theories [20T4221 ] , and it is how the initial curvature 
perturbations are created in the curvaton model [2314261 ] . 
In this case, the second scalar field's decay supersedes 
the inflationary reheating epoch, and its fluctuations set 
the spectrum of adiabatic perturbations. Clearly, this 
process could occur multiple times, with the caveat that 
thermal relics may survive subsequent reheatings. 

The final round of reheating takes particular impor- 
tance because it generates the thermal plasma of the Uni- 
verse and fixes the spectrum of perturbations on scales 
greater than the horizon size during the last reheating 
epoch. Physical processes that occur at temperatures 
larger than the temperature of the radiation bath fol- 
lowing the last reheating epoch (the reheat temperature) 
will be affected by the altered expansion history of the 
Universe prior to reheating. The impact of a low re- 
heat temperature on the relic abundance of WIMP dark 
matter has been studied extensively [27l - l33 | ; a low reheat 
temperature reduces the relic abundance of thermal dark 
matter, but the abundance may be enhanced if the dark 
matter is also produced nonthermally. A low reheat tem- 
perature also relaxes constraints on axion [2l| [3|1 [36| and 
neutrino [29l.l37l.l38j dark matter. This prior work on low- 
temperature reheating has only considered its impact on 
the homogeneous content of the Universe, but deviations 
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from radiation-dominated expansion also affect density 
perturbations. In this paper we consider how the evo- 
lution of the Universe prior to the creation of the final 
radiation bath influences perturbations on scales that en- 
ter the horizon before the last reheating epoch. 

A homogeneous scalar field that rapidly oscillates 
around the minimum of a quadratic potential has the 
same dynamics as a pressureless fluid of non-relativistic 
matter — its energy density scales as a -3 |l4|. More- 
over, it has been shown that the correspondence between 
a rapidly oscillating scalar field and a matter fluid ex- 
tends to the perturbations (39l - l4^ |. In particular, sub- 
horizon perturbations in the scalar field's energy density 
with fcphy S < y^SHm^, where is the mass of the scalar 
field, grow linearly with the scale factor while the scalar 
field dominates the Universe [4l|, H^]. If the scalar os- 
cillates for a sufficiently long time before decaying, den- 
sity perturbations in the scalar field become nonlinear 
and could produce gravitational waves [HI, Hi[ . Longer 
wavelength perturbations enter the horizon later and re- 
main linear prior to reheating. Nevertheless, the growth 
of these perturbations during an early scalar-dominated 
phase may have a significant impact on the subsequent 
growth of structure. 

To assess effects of an early "matter" -dominated epoch 
on the present-day matter power spectrum, we study 
the evolution of perturbations through the reheating era. 
Our analysis applies to any scenario in which the energy 
that dominates the Universe prior to reheating behaves 
as a pressureless fluid and decays perturbatively. This 
fluid could be composed of metastable non-relativisitic 
particles, or it could be an oscillating scalar field. We 
note that a coherently oscillating scalar field may decay 
non-perturbatively; during preheating, the oscillations of 
the scalar field lead to resonant particle production that 
cannot be captured by standard perturbation theory (TB - 
fl7j | . A coherently oscillating scalar field still decays per- 
turbatively, however, if its couplings to other fields are 
sufficiently small that preheating is unimportant. Specifi- 
cally, our analysis applies to the decay of coherent scalar 
fields in the narrow resonance regime; due to the Uni- 
verse's expansion, preheating is inefficient in this regime, 
and the scalar 's decay can be treated perturbatively [l7| . 

Inspired by models in which dark matter is a d ecay 
product of a string modulus [H, HI], the inflaton (UHlf, 
or the curvaton [49( , we allow the pressureless fluid in our 
reheating model to decay into cold dark matter particles 
and radiation. We focus on modes that are linear during 
scalar domination; these are the largest scales that enter 
the horizon prior to radiation domination and thus have 
the highest potential for observational impact. We find 
that the radiation perturbation retains no memory of the 
scalar perturbation's growth. On the contrary, the radi- 
ation density perturbation on scales that enter the hori- 
zon during the "matter" -dominated phase are dramat- 
ically suppressed compared to perturbations that enter 
the horizon after reheating. Consequently, dark matter 
that is coupled to or produced from the radiation after 



reheating will have a cutoff in the matter power spectrum 
inherited from the last reheating epoch. Alternatively, if 
nearly-decoupled dark matter is produced directly at re- 
heating, the dark matter inherits the enhanced density 
perturbation that grew prior to reheating. In this case, 
we get an enhancement of the small-scale matter power 
spectrum and an epoch of early small-scale structure for- 
mation during the cosmic dark ages. 

Our paper is organized as follows: In Sec.|TI]we discuss 
our three-fluid model for reheating and its background 
evolution, and we present the equations that govern the 
evolution of the fluid perturbations and their initial con- 
ditions. In Sec. Mil we describe the evolution of the per- 
turbations in the radiation fluid through reheating. In 
Sec. II VI we describe the evolution of the cold dark mat- 
ter decay product, including the effects of free-streaming, 
and we derive a new transfer function for the matter per- 
turbation. In Sec. [V] we show that the new small-scale 
matter power spectrum leads to an abundance of low- 
mass structures at high redshift. In Sec lVII we summarize 
our results and briefly discuss the prospects for detecting 
the numerous compact dark matter subhalos predicted 
by this reheating scenario. Appendix [X] contains the 
derivation of the perturbation evolution equations, and 
the initial conditions for the perturbations are derived 
in Appendix [B] Finally, Appendix [C] provides expres- 
sions necessary to evaluate the matter transfer function, 
including a scale-dependent growth function for small- 
scale perturbations. 



II. THREE-FLUID MODEL FOR REHEATING 

Our reheating model consists of a pressureless fluid 
that dominates the energy density of the Universe be- 
fore it decays into radiation and dark matter particles. 
We assume that the radiation is tightly coupled and be- 
haves like a perfect fluid with P x = p r /3. We also as- 
sume that the dark matter particles are nonrelativistic 
and completely decoupled from the radiation bath from 
the moment of their creation. We initially neglect the ve- 
locity dispersion of these particles, but in Section IIVBI 
we consider the impact of free-streaming by dark matter 
particles that are created with a nonzero velocity. 

Inspired by the curvaton and inflaton, and to avoid 
confusion with dark matter produced during reheating, 
we call the pressureless fluid that initially dominates the 
Universe a scalar field. As mentioned in the previous sec- 
tion, an oscillating scalar field behaves like a pressureless 
fluid, with p oc a~ 3 , if the period of the scalar oscilla- 
tions is much shorter than the Hubble time (m^ ^> H). 
We will assume that the decay of this scalar field is slow 
enough that the coherence of the scalar oscillations are 
irrelevant; in the language of preheating, we are assum- 
ing that the scalar decay is the narrow resonance regime 
in an expanding universe. If the scalar decay rate 
satisfies T^/m^ <C (m^/mpi) 2 , then momentum modes 
of the scalar's decay products redshift through their in- 
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stability bands too quickly to be significantly excited. 1 
In this case, the decay of the scalar field can be treated 
perturbatively and the scalar's coherence is irrelevant. 

In this three-fluid model for reheating, the background 
equations for the energy density p$ of the scalar field, the 
radiation density p r , and the dark matter density pdm are 



dt 



dt' 
d_ 
dt 

Pdn 



+ 3Hp<p = -T^p,),, (1) 

p r +AHp r = +(1-/)I> , (2) 

3Hp dm = d-fY^p^, (3) 



where / is the fraction of the scalar's energy that is trans- 
ferred to dark matter particles. For a given reheat tem- 
perature, there is one value of / that gives the observed 
dark matter abundance, as discussed below. In these 
equations, we have considered only interactions that are 
linear in the scalar field; such interactions are required 
for the scalar to decay completely fl8j . Nonlinear in- 
teractions, which correspond to scalar self-annihilations, 
would introduce terms proportional to p^ in the equa- 
tions above. If one considers an oscillating scalar field 
that decays into a massless scalar field \ y ia an interac- 
tion term proportional to </>x 2 , the Boltzmann equation 
averaged over many oscillations implies that p x follows 
Eq. © [gjj. 

Similar three-fluid models for reheating have been used 
to calculate the relic abundance of dark matter in low- 
reheating-temperature cosmologies [HI, [H, IHj]. These 
models also include dark matter self-annihilations and 
thermal production. By neglecting such interactions in 
Eqs. © and ([3]), we are assuming that the dark matter 
particles produced in scalar decays far outnumber any 
dark matter particles that may have been produced ther- 
mally. We are also assuming that the velocity-averaged 
self-annihilation cross section ((av)) of the dark matter 
particles is small enough that the nonthermal dark mat- 
ter does not self-annihilate. For any value of (av), these 
assumptions are valid for sufficiently low reheat tempera- 
tures. If (av) is too large, however, these reheat temper- 
atures may be below 3 MeV and are therefore disallowed 
by BBN 

With these assumptions, the duration of the radiation- 
dominated era after reheating is determined by the 
branching ratio / and does not depend on (av) [32l |. 
If ?rh is the reheat temperature, and T cq is the 
temperature at matter-radiation equality, then / ~ 
0.43(T cq /T RH )(10.75/ 5 * s ) 1/3 , where g, s is the entropy 
density divided by (27r 2 /45)T 3 and should be evalu- 
ated at the reheat temperature. 2 From the cosmic mi- 



1 This condition follows from demanding that q >C H/m^ <C 1, as 
in Ref. [13 

2 We focus on non-relativistic dark matter at reheating. Models 
with relativistic dark matter at reheating [5l| need a higher / for 
the same late-time abundance, but we do not consider this case. 




scale factor (a) 

FIG. 1: The evolution of the energy densities in the decay- 
ing scalar field and its decay products, as fractions of p C rit,o: 
the critical density at a = 1. The scalar field is modeled as 
a pressureless fluid; its energy evolves as p cc a -3 until the 
expansion rate equals its decay rate at a ~ 10 4 . While the 
Universe is scalar dominated, the scalar is feeding energy into 
the matter and radiation fields, and both the matter and ra- 
diation fields evolve as p oc a~ 3//2 . Once nearly all the scalar 
density is transferred to the matter and radiation fields, their 
energy densities evolve in the usual way: radiation as p cc a~ 4 
and matter as p oc a -3 . 



crowave background, we know that T cq = 0.75 eV [10(, so 

/ ~ (6.4 x Ifr 8 )(5 MeV/T RH )(10.75/ 5 * 5 ) 1/3 is required 
to produce the observed amount of dark matter. 

When numerically solving Eqs. CO)-© for the evolu- 
tion of these three fluids and their perturbations, we de- 
fine an initial time to with a(irj) = 1 and H(to) = H±. At 
t = to, the critical density is pcrit.o = 3mp[i/ 2 /(87r), and 
we define dimensionless density variables p^ = p^/ p cr itfl, 
Pdm = Pdm/Pcrit,o, and p r = Pr/Pcxitfi- We also define 
a dimensionless decay rate = Y^/H\. We will as- 
sume that the scalar field initially dominates the energy 
density of the Universe. The scalar field decays when 
H ~ r^,. Since the Universe is effectively matter domi- 
nated prior to the decay of the scalar field, H oc a~ 3 / 2 , 
and the value of the scale factor at the time of scalar 

~ 2 /3 

decay is a R u — .At this time, the Universe transi- 

tions from scalar domination to radiation domination. 

In the limit that Y^t -c 1, the Universe is scalar dom- 
inated, and the solutions to the background equations 
are 



P4> = P4>(t )a~ 



Pdn 



= ^^Mto)iy ]a- 3/2 

= /MMlWo]^ 372 +Cr, 



+ c r a 

-3 



(4) 
(5) 
(0) 



In these equations, c r and c m are arbitrary constants to 
be determined by initial conditions. If there is radiation 
or dark matter that does not originate from scalar decay, 
then c r and c m are positive. If there was a time when 



4 



there was no radiation or dark matter, then c r or c m are 
negative to make p r = or pd m = at that time, as in 
Ref. (29[. In either case, the contributions to the radia- 
tion and matter densities from the terms proportional to 
c r and c m become less significant as the scalar-dominated 
era continues. We set our initial condition late in the 
scalar-dominated era, after any matter or radiation that 
did not originate from scalar decay is diluted and the 
system loses its memory of the beginning of scalar de- 
cay. In this case, c r = c m = 0, and the initial values 
(at t — to) of p~0,p r , and pdm are determined by the 
chosen values for and the branching ratio /. Dur- 
ing scalar domination, H = 2/ (3£), which implies that 
P0(£o)IWo = (2/3)r <j> p^{t ). It follows that 

2 - 

Pdm (^o) = g./T0P0(to), (7) 
Pr(to) = jj(l-/)f p (t o ). (8) 

Finally, our assumption that the Universe is flat demands 
that P(j>{to) + p r (to) + Pdm{to) = 1. We want the Universe 
to be initially dominated by the scalar energy density, so 
we must set -C 1. Figure Q] shows the evolution of the 
energy densities of the scalar, matter and radiation fields 
for f = 1(T 6 and / = 1(T 8 . We see that the Universe 

~ 2/3 

becomes radiation dominated at orh — Fj . 

As described in Appendix [XJ the evolution equations 
for perturbations in the scalar, radiation and matter flu- 
ids are derived by perturbing covariant versions of Eqs. 
(HJ-© [55 - [54j . We work in conformal Newtonian gauge, 
with 

ds 2 = a 2 (r)[-(l + 2*)dr 2 + %(1 + 2$)dx i dx j ]. (9) 

In addition to the metric perturbations 9 and $, 
we must solve for the fractional density perturbations 
[8 = Sp(r, x)/p°(t)] in the scalar field (8<p), the radia- 
tion (S r ), and the matter (Sd m )- The three fluids also 
have velocity perturbations [v* = dx l /dr], and we solve 
for the divergences of the velocity fields for the scalar 
field (00 = V ■ Vcj,), the matter (9d m = V • ildm) and the 
radiation (0 r = V • v r ). The perturbation equations are 

8$ + + 3$ = -aT^, (10a) 
00 + ^0 + V 2 * = O, (10b) 
• 4 • pl 

S r + -6 r + 4$ = (1 - f)aY^ [8$ - S r + 9] , (10c) 
3 p r 

Or + V 2 (| + *) = (1 - f)aT/^ - 0.) , (lOd) 

P° 

Sdm + 0d m + 3$ = JaT^ % - <5 dm + ¥] , (lOe) 

Pdm 



6dm + V 2 * + -0dm = faT^ [00 - 9dm] , (10f) 
a Pdm 



where a dot denotes differentiation with respect to con- 
formal time t. We also have the perturbed Einstein field 
equation 

V 2 $+3^ (^B - $j = -AnGa 2 (p ^ + p% + Pd m <W) • 

(11) 

Finally, we close the system of perturbation equations 
by noting that the absence of anisotropic stress implies 

To numerically solve these equations, we go to Fourier 
space and define the following dimensionless parameters: 
E(a) = H(a)/H lt k = k/H u 00 = 9^/ Hi, 9 r = 9 r /H lt 
and 0dm = 0dm/ H±. The equation set that we solve is 
given in Appendix [B] We use the scale factor a as our 
time variable, and we set initial conditions when the Uni- 
verse is scalar dominated and the mode is outside the 
Hubble horizon. Since we start our solution when a = 1 
and H = Hi, these restrictions imply that r0 -C 1 and 
k -C 1. In Appendix [Bl we derive the following initial 
conditions: 



$(a ) 


= $0 


(12a) 


8<t>{ao) 


2~9 

= 2$ + gAr$ Q Oo 


(12b) 


8 r (a ) 


= $a + —k z <f> a Q 
63 


(12c) 


<W(ao) 


2~9 

= $0 + gAr$ ao 


(12d) 


(ao) = 0r(a o ) 


2 _ 

= 0dm(a o ) = --fc 2 $v^o. 


(12e) 



with ao = 1. These initial conditions neglect terms that 
are O(k 4 a ) and 0(r0ao). 

The initial conditions presented above differ from the 
characteristics of superhorizon adiabatic perturbations in 
a radiation-dominated universe, for which 8 r = 2$ and 
8dm = (3/4)5 r . Since the matter and radiation origi- 
nated from a single scalar field, we should be left with 
adiabatic perturbations in the matter and radiation af- 
ter the scalar vanishes. Indeed, modes that are super- 
horizon at the time of reheating [k < Qrh-£' (^rh)] evolve 
when a ~ cirh, and this evolution takes the perturbations 
from the initial values given above to the conditions for 
adiabatic superhorizon modes in a radiation-dominated 
universe. The gravitational potential transitions from its 
initial value $0 = $(*o) to (10/9)<&o as the Universe tran- 
sitions from being effectively matter dominated to being 
radiation dominated. The radiation density perturbation 
grows from S r (t ) = $0 to 2(10/9)$ O - Meanwhile, the 
matter density perturbation grows from Sdm(to) = $0 
to (5/3)$o- Thus, after the Universe becomes radia- 
tion dominated, we have <5dm = (3/4)<5 r on superhori- 
zon scales. Once these modes enter the horizon, they 
behave exactly as expected for adiabatic perturbations 
in a radiation-dominated universe, with $ = (10/9)<i>o 
on superhorizon scales. In the following sections, we will 
consider the smaller-scale modes that enter the horizon 
while the Universe is still scalar dominated. 
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III. PERTURBATIONS IN THE RADIATION 
FLUID 

We briefly review the standard evolution of a perfect- 
fluid radiation perturbation in a radiation-dominated 
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scale factor (a) 

FIG. 2: Evolution of <5 r /$o , &r/&o, and $/$o for a mode 
that enters the horizon at a ~ 1000. In this cosmological 
scenario, the Universe is radiation dominated; there is no era 
of scalar domination. 



where $o is the initial value for $ set when the mode is far 
outside the horizon (55[, and x = kr — ka/[a^ t:D E(anx))]. 
Here, ord denotes a value of the scale factor well after 
reheating (deep into radiation domination) and ^(ord) 
is the dimensionless Hubble factor as discussed in Sec. HT1 
In Fig. [2] we show the evolution of the radiation pertur- 
bations in the standard scenario. The perturbations are 
frozen when a given mode is superhorizon and oscillate 
after the mode enters the horizon. As the amplitude of 
oscillations in $ decay as A$ oc x~ 2 ■> the equations for 
the perfect-fluid variables S r and 9 r reduce to those of 
a harmonic oscillator. After they are well in the hori- 
zon the perturbations in S r and 9 r quickly asymptote to 
cosine and (out-of-phase) sine solutions with amplitudes 
Ag r = 6$o and Ag = (k-^3/4)As r respectively. 

We now compare the subhorizon behavior of the ra- 
diation perturbations following a period of reheating to 



universe, and then compare this to the evolution of radi- 
ation perturbations when followed through reheating. 

In a radiation-dominated universe, the equations de- 
scribing the evolution of perfect-fluid radiation pertur- 
bations can be solved exactly and take the form 



(13a) 
(13b) 
(13c) 

I 

the standard behavior given in Eq. (|13[) . To make this 
comparison, we must take care to compare modes that 
have the same physical wave vector at some temperature 
Trd long after the reheating process is complete. These 
modes have not typically spent the same amount of time 
inside the Hubble horizon. If a mode enters the horizon 
prior to reheating, then it does so earlier than it would 
have if the Universe never had an epoch of early scalar 
domination. For instance, the mode shown in Fig. Ren- 
ters the horizon at a = 1000 if the Universe was always 
radiation dominated. But if the Universe was scalar dom- 
inated up to a ~ 10 , then this mode enters the horizon 
when a = 100. In general, for modes that enter the hori- 
zon during scalar domination, the value of a at horizon 
entry in the scalar-dominated scenario is fcRn/fc times the 
value it would have been in the standard scenario, where 
fcRH = clrhH (orh) is the wave number of the mode that 
enters the horizon at reheating. 

Figure [3] shows the same mode as Fig. [5] in a universe 
that is scalar dominated for a <S 10 4 (corresponding to 
the choice = 10~ 6 ). The evolution of this mode is 
drastically different than the standard evolution shown 
in Fig. [3J Initially, 5 r grows considerably while the Uni- 
verse is scalar dominated. However, after reheating this 
enhancement is lost and 6 r oscillates with a small am- 
plitude. The late-time oscillations have the same pe- 
riod as the oscillations given by Eq. (fT3"|). but the am- 
plitude of the fluid perturbations is suppressed by a fac- 
tor of ~70.6 with respect to the standard case, with 
A Sr ~ (6/70.6)$o - 0.085$ and A §r = {kV3/A)A Sr . 
The $ solution does not begin to decay until after reheat- 
ing in this scenario but subsequently decays as x~ 2 with 
a similarly suppressed amplitude relative to the standard 
evolution. Figure U shows the k = 0.1 mode in a universe 
with = 10~ 9 . The late-time solution exhibits simi- 
lar behavior, only now with an amplitude suppressed by 
^7600 with respect to the standard case, so that A$ r ~ 
(6/7600)$o - 0.0007$o- In short, the radiation pertur- 
bation at late times follows, up to a phase shift, Eq. (fl~3"|) 



9£o 

X 3 

X 3 



VSsin ( —= ) 
n/37 



X cos 



\V3 



2V3( X 2 - 3) sin ( JQ - X (X 2 - 6) cos 



2v / 3xcos 



(x 2 



6)sin (^) 



6 





V*o 






VfH^o) 




- 













10° 10 1 10 2 10 3 10 4 10 5 



scale factor (a) 

FIG. 3: Evolution of 5 r /$o, r /$o, and $/$o for a mode 
(k — 0.1) that enters the horizon at a ~ 100. In this cosmolog- 
ical scenario, the Universe was scalar dominated for a ^ 10 4 , 
and then it became radiation dominated. When the Uni- 
verse is radiation dominated, <5 r oscillates with an amplitude 
of 0.085$o- For this mode, fc/fcRH — 11. 
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FIG. 4: Evolution of <5 r /$o, 6 r /<&o, and <E>/<E>o for a mode 
(fe = 0.1) that enters the horizon at a ~ 100. In this cosmolog- 
ical scenario, the Universe was scalar dominated for a ^ 10 6 , 
and then it became radiation dominated. When the Uni- 
verse is radiation dominated, 5 r oscillates with an amplitude 
of 0.0007$ . For this mode, fc/fc RH ^ 114. 



with a suppressed amplitude that depends on the value 
of fc/fcRH- For fc/fc RH £ 20 we find A Sr < 10" 3 $ with 
the precise value modulating by an order of magnitude as 
a function of fc/fcRH- For 20 <; fc/fcRH *t 3 A,5 r increases 
with decreasing fc, reaching a peak of A,^ ~ 10$o- For 
smaller values of fc, A$ r decreases as k decreases and flat- 
tens out at a value of A$ r ~ (60/9)$o f° r fc/fcRH ^ 0.1. 
Modes with fc/fcRH < 0. 1 are superhorizon at the time of 
reheating and thus their evolution follows Eq. (fl"3| with 
$0 — > (10/9)$o since in this scenario $o is defined during 
a period of matter domination. 



This behavior can be understood by re-examining 
Eqs. (|10c[) and (|10dj) . During the scalar-dominated 
epoch, $ = 0, 8^ cx a, and 5$ grows large. Since the 
scalar is decaying to radiation, this large scalar over- 
density becomes a source for 8 r , and 8 r grows consid- 
erably. The increase in the radiation overdensity 8 r in 
turn sources a higher 9 r , representing an outflow of radi- 
ation fluid in real space, through Eq. (|10dl) . This outflow 
slows the growth in S r , which eventually saturates and 
has 8 r = 0. Meanwhile, as the decay continues, 9 r keeps 
growing and the outflow balances the source terms from 
scalar decay in Eq. (llOcj) . Eventually, the scalar source 
term becomes subdominant, and 8 r decreases as the out- 
flow continues. This decrease in 8 r slows the outflow, and 
9 r decreases to a point where 8 r and 8 r begin to oscil- 
late out of phase. Due to the relatively large value of 
8$, these oscillations are forced; the equilibrium point is 
shifted to positive 5 r , and the first trough in S r is very 
shallow. However, by the time S r emerges from the first 
underdensity, the scalar density is zero and henceforth 8 r 
undergoes free oscillations with a small amplitude. 

In summary, the radiation perturbation long after re- 
heating is described by Eq. |[T5]) , but with <J> replaced by 
T(fc)$o and a — > a + arj, where arj is the phase shift re- 
sulting from the fact that oscillations started at reheating 
rather than horizon entry. For fc <; 20fcRH, T(k) ^ 10~ 3 , 
and the perturbations are suppressed. As fc decreases, 
T(fc) increases, reaching T(fc) ~ 1.5 for 2 < fc/fcRH ^5 4. 
For smaller values of fc, T(k) again decreases until it lev- 
els out at T(fc/fc RH & 0.1) = 10/9. The major effect of 
a period of scalar domination on the radiation perturba- 
tion is a suppression of power on scales that entered the 
horizon during the scalar-dominated era. These scales 
are very small as /crh is given by 




where Trh is the temperature of the radiation bath when 
the Universe becomes radiation dominated. In this ex- 
pression, 5* = p r /[(7r 2 /30)T 4 ], and both and g„s are 
evaluated at Trh- If the dark matter is a thermal relic 
that decoupled from the radiation bath after reheating, 
perturbations with fc > /crh will be suppressed. 

This suppression of density perturbations is in prin- 
ciple relevant for any theory in which dark matter has 
origins in the thermal bath after reheating, such as the 
standard WIMP paradigm. The characteristic kinetic 
decoupling temperature Trd in WIMP models is nec- 
essarily below the chemical decoupling temperature and 
thus the reheat temperature Trh if they are thermally 
produced after reheating (56|. However, in models where 
kinetic decoupling occurs immediately after chemical de- 
coupling, the nominal cutoff from kinetic decoupling and 
the cutoff from altered perturbation evolution in this re- 
heating scenario can be comparable and a more detailed 
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FIG. 5: The evolution of the fractional density perturba- 
tion in the matter field (5dm)- This mode enters the horizon 
when a = 25, which is during the scalar-dominated era. Dur- 
ing scalar-domination, subhorizon matter perturbations grow 
linearly with a. The Universe becomes radiation dominated 
when a = 10 4 , at which point 5dm grows logarithmically. The 
long-dashed curve depicts Eq. (|15p . which describes the evo- 
lution of (5dm during scalar domination. The short-dashed 
curve depicts Eq. (|16|) with qrh = 1.29f^ 2//3 , which fits the 
evolution of 5dm after reheating. 



analysis including both effects is warranted. In mod- 
els where the dark matter is never in equilibrium with 
the standard model but is nevertheless produced from a 
thermal plasma after reheating, such as in feebly inter- 
acting massive particle (FIMP) models where the dark 
matter abundance freezes- in [57jj , in h ylog enic models 
[l2} . in some hidden-sector models [HI, l58j ] and gener- 
ally for models where Tkd ^ Irh, the reheating cutoff 
is expected to be the dominant cutoff in the small-scale 
matter power spectrum. 



<5dm- Prior to the decay of the scalar field, the Universe 
is effectively matter dominated; a 2 E(a) — y/a and $ re- 
mains constant after the mode enters the Hubble horizon. 
The solutions derived in Appendix [B] for the velocities 
9,p(a) and #dm(a) hold until the Universe becomes radia- 
tion dominated. It follows that 



S 'dm( a ) 



;fc 2 $o 



3 

2a 



-Mm 



$] . (14) 



The quantity in the square brackets is initially zero, 
and it remains zero until the Universe becomes radiation 
dominated, as one can verify by considering S^ — (5 dm . 
Therefore, during scalar domination, 



2 - 

Sdm(a) = $0 + gfc 2 $oa, 



(15) 



and we see that the early-time solution derived in Ap- 
pendix[B]holds after the mode enters the Hubble horizon. 
This model for <5 dm (a) is depicted in Fig. [5J and we see 
that it matches the numerical solution for a ^ csrh- 

The linear growth begins when the linear term ex- 
ceeds the initial value [(2/3)fc 2 a <; 1]. We define 
aijn = 3/(2fc 2 ). For modes that enter the hori- 
zon during scalar domination, cthor = k , so ay ln 
is larger than cihor- Since Sdm = $0 prior to the 
onset of linear growth, at the time of scalar decay 
<W(grh) = (aRH/aiin)$o = (2/3)fc 2 a R ,H"I > o- After the 
Universe becomes radiation dominated, (5dm (a) grows log- 
arithmically, and a<5 dm (a) is constant. We can therefore 

set a5' dm (a) = aR H <5 dm (aR H ) = (2/3)fc 2 i> aRH- It follows 
that, for modes with k/kjm <; 1, 



<5dm(a) = -a RH fc 2 <E>o 



1 



In (-!- 

\orh 



(16) 



after reheating. This model for (5 dm (a) is depicted in 
Fig. [5J and we see that it matches the numerical solution 
for a £ cirh. 



IV. PERTURBATIONS IN THE MATTER 
FLUID 

Since the radiation perturbation does not retain any of 
its growth during the scalar-dominated era, an epoch of 
scalar domination results in primordial structure growth 
only if the dark matter decouples from the radiation be- 
fore the Universe becomes radiation dominated. We now 
consider the evolution of the decay-produced matter per- 
turbation Sdm for modes that enter the horizon while the 
Universe is dominated by the oscillating scalar field. For 
these modes, (5 dm grows linearly with a between horizon 
entry and reheating, and then Sdm grows logarithmically 
after the Universe becomes radiation dominated. An ex- 
ample is shown in Fig. [51 which shows the evolution of 
<5dm for k = 0.2 and = 10~ 6 ; this mode enters the 
horizon when a = 25, and the Universe becomes radia- 
tion dominated when a = 10 4 . 

We now derive an analytic model for the evolution of 



A. The Matter Transfer Function 

Due to the superhorizon evolution of $ during the tran- 
sition from scalar domination to radiation domination, 
the definition of the transfer function is ambiguous. It 
is standard to define the transfer function T{k) through 
the relation 

3 k 2 

<5 dm (fc, a > a oq ) - -t^—fti ^ P {k)T(k)D(a), (17) 

where Hq is the present-day Hubble parameter, Qm is the 
current density of matter divided by the current critical 
density, <f> p is the initial potential fluctuation, and D(a) 
is the growth function normalized to make D(a) = a dur- 
ing matter domination. With this definition, T{k) — 1 on 
large scales because modes that enter the horizon after 
matter-radiation equality have $ = (9/10) $p at hori- 
zon entry. In our scenario, however, these modes enter 
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FIG. 6: The fractional perturbation in the matter density 
evaluated after the scalar decays (a = lOOOf 7 2//3 ). The solid 
curve is the numerical evaluation. Modes with k/k^M < 0.001 
are still superhorizon at this value of the scale factor. The 
long-dashed curve is a plot of Eq. (|18p and describes subhori- 
zon modes that enter the horizon after the scalar decays. The 
short-dashed line is a plot of Eq. (|19p and describes modes 
that enter the horizon during scalar domination. 



Modes with fc/fcRH <i 1 enter the horizon while the 
Universe is dominated by the oscillating scalar field. To 
obtain the transfer function for these modes, we express 
<5dm(a) in the same form as Eq. (|18p and then extract 
the values for A and B. From Eq. (|16l) . we see that 
A = (3/5)aRnfc 2 and ln£? = 1 + ln(ahor/aRH)- We now 
evaluate A and B in terms of k/k^ji- First we note that 
we never used orh = f ^ 2 ^ 3 in the derivation of Eq. (|16l) . 
Rather, orh was used as the transition point between 
linear and logarithmic growth. We find that Eq. (fT!)]) fits 

the numerical solution for <5dm(a) if Orh = l^f^ 2 ^ 3 , 
as shown in Fig. [5] We also find that f V 3 ~ 1.29fc RH 
for a wide range of values of T^. We conclude that we 
should use o,rh = k^ when evaluating A and B. Since 
these modes enter the horizon deep in the era of scalar 
domination, we also have ahor = k~ 2 . We conclude that 



Sdm (a) 
A 
B 



10 
~9 

3 
5 



A In 



Ba 

G>\iov 



(19) 



k 

&RH 



l RH 



the horizon with $ = $o> where $o is the superhori- 
zon potential perturbation during scalar domination. To 
define our transfer function, we will use the standard 
definition above, but we note that <£> p is the superhori- 
zon potential perturbation during radiation domination; 

$ p = (10/9)$„. 

Figure [6] shows how Sdm at a set value of the scale 
factor depends on k. In the cosmology depicted in this 
figure, = 10~ 6 , which implies that <2rh = 10 4 and 
A;rh = aRH^(flRH) = 0.0088. The matter perturbation 
is numerically evaluated at a = 10 7 for a wide range 
of k values; this is the solid curve in Fig. [5] We see 
that the function (5d m (1000aRH, k/k^m) contains three 
different behaviors. First, for fc/fcRH ^ 0.001, we have 
<5dm(1000aRH, k/ksji) — 5/3$o- These modes have not 
yet entered the horizon at a = 10 7 , and so they are fixed 
at the adiabatic "initial" condition, as described in Sec- 
tion nu 

Modes with 0.001 ^ fc/fcRH & 1 enter the horizon dur- 
ing radiation domination. These modes grow logarithmi- 
cally once they enter the horizon; as long as Sp r ^ Spdm, 
their behavior is well-described by the function 



<5dm(a) = — $0 

y 



A In 



Ba 

Qhor 



(18) 



with A = 9.11 and B = 0.594, as given in Ref. [59j . 
The long-dashed curve in Fig. [5] shows this function with 
a = 1000<2rh and a\ lor E(a\ lo - l: ) = k. We see that it is an 
excellent fit for 0.01 ^ k/kmi & 1- We conclude that the 
standard transfer function can be used for modes with 
fc/fcRH 1; these modes are unaffected by the era of 
scalar domination. 



for modes with fc/fcRH <; 1 evaluated after scalar de- 
cay while Sp r ^> Spdm- This model, with a = 10 7 and 
ahor^(ahor) = k, is shown by the short-dashed line in 
Fig. |6] and we see that it is an excellent fit to the nu- 
merical solution when fc/fcRH <; 10. To smoothly con- 
nect this solution to the k ^ /crh solution (A = 9.11 and 
B = 0.594), we found functions A{k/k^s) and B{k/ks.s) 
that fit the numerical solution for <5dm(a, k) when they 
are inserted into Eq. (|19j) . These functions are given in 
Appendix [Cj 

Now that we have A(k) and B(k), we can obtain the 
behavior of (Sdm (a) during the matter-dominated era by 
matching Eq. (|T9]) to the decaying and growing modes of 
the Meszaros equation, which is valid when Spdm ^ Sp r 
60]. While the baryons are still coupled to the photons, 
they do not fall into the potential wells created by the 
dark matter density perturbations. Prior to baryon de- 
coupling and after matter-radiation equality, Sd m (a) is 
given by !59] 



. 3A {10 

Sdm = —fl I ln 



4\ft Ba 



eq 



Oho 



V (a) , (20) 



where f± , ji and a are determined by the baryon fraction 
fb = Pb/(Pb + Pdm), 

h = 1 - 0.568/ b + 0.094/ 2 

h = l-1.156/ b + 0.149/ b 2 - 0.074/^, 

and T>(a) is the growing solution to the Meszaros equa- 
tion. If baryons do not participate in gravitational col- 
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lapse, 







— a 






a Gq J 


2^1 


a, a 



1 a, 



cq 



2 a + a, 



cq 



(21) 

where 2-F1 [a, b; c; ac] is Gauss's hypergeometric function, 
and 



a= 4 



1 - VT + 24(1 - /b 



(22) 



Long after matter-radiation equality, Eq. ([2~Tj) reduces 
to 2?(a) = (a/a oq ) _Q ; dark matter overdensities experi- 
ence slower than linear growth if the baryons do not fall 
into their potential wells. Consequently, matter pertur- 
bations on scales smaller than the baryon Jeans length 
grow slower than larger-scale perturbations [55j , and we 
will need to account for this suppression in our transfer 
function. 

To complete our derivation of the matter transfer func- 
tion, we evaluate a eq /ahor for modes that enter the hori- 
zon during scalar domination: 



y/2k ( 



icq 



9*s[T] 



RHJ 



3.91 



2/3 



3.36 



5*Pi 



RH 



1/2 



(23) 

The terms containing and g*s account for changes in 
the number of relativistic species between reheating and 
matter-radiation equality, and the product of these terms 
is always between 1.1 and 1.6. The product Ba eq /ah T is 
proportional to kmi/k eq , and it appears in Eq. (|20[) log- 
arithmically. Since fcptH/fccq 3> 1, the terms containing 
5» and g*s have a negligible effect on Sdm, and we do not 
consider them further. For modes that enter the horizon 



during radiation domination, 
in general, we find that 



j/ahor = V2fc/fc cq , and 



^cq 
^hor 



V2k 



■"■cq 



1 ( k V 

V &rh J 



4.2351 1 / 4 - 235 



(24) 



for all 



fits the numerical solution to k = ah m H(cih, 
values of k. 

Inserting Eq. (|23|) into Eq. (|20|) and evaluating the 
transfer function without baryons gives 



T(k » A; RH ) 



RH 



In 



eV2k RH 

P 3 



k, 



cq 



(25) 



Thus we see that T(k ^> /crh) is scale-invariant. These 
modes grow linearly with the scale factor between hori- 
zon entry and scalar decay, and they grow linearly 
with the scalar factor after the Universe becomes mat- 
ter dominated. The transfer function characterizes de- 
viations from linear growth in 5a m (a). During radia- 
tion domination, these modes grew logarithmically in- 
stead of linearly, and so the transfer function depends 
only on the duration of the radiation-dominated era: 



T(k > fc RH ) ~ (k-RH/k C q) 2 ln[fc RH /fc eq ]. Given current 
measurements (k eq = 0.0098 Mpc~ 4 ) we note that 



^ = 1.72 x 10 11 ( Trh 



VlOOGeV; V.9*5 / V100/ 



1.18 x 10 b 



RH 



lMcV 



10.75\ 1/3 f g* 



9*s 



10.75 



1/2 



where g* and g*s are evaluated at Trh. 

When we consider structure formation in SectionfVl we 
will need a transfer function that is applicable to all k val- 
ues, including k <5 fc eq , and includes the effects of baryons 
and neutrinos on structure formation. We use CAMB 
Sources [6l| to compute the matter transfer function for 
k/k cq < 8.2 x 10 5 . On very small scales {k/k eq <; 5x 10 5 ), 
the transfer function computed by CAMB (Tcamb) has 
the same scale dependence as the transfer function (Teh) 
provided by Eisenstein & Hu (1998) [62], so we extend 
the transfer function to larger k by taking 

Wfc/fc ^BOvlf^ T ,,^ T CAMB(fc/fc cq ^8.2x 10 5 ) 

T (k/k cq > 8.2 x 10 ) = TMk) TMk/kcq = 8 _ 2xlQ5) • 

When evaluating the transfer function, we take 
H = 70.4 km/s/Mpc, Q M = 0.272, and fl h = 0.0456, 
following Ref. [Ia |. 

Even after recombination, the baryons have nonzero 
pressure, and they do not participate in gravitational col- 
lapse on scales that are smaller than the baryon Jeans 
length 1/fcjb [H|. While matter perturbations on scales 
with k < fcjb grow linearly with the scale factor after 
recombination, the growth of smaller-scale perturbations 
still follows T>(a), given by Eq. (|2~T1) . We account for 
this suppression by using CAMB Sources to determine 
the density perturbation at a redshift of 50, and then 
we use Eq. (fTT|) to find T(k)D(z = 50), which we re- 
fer to as the transfer function at z = 50. We then 
apply a scale-dependent growth function, normalized to 
unity at z = 50, to obtain the density perturbation at 
other redshifts. At z = 50, Tcamb (k) — Teh(^) for 
k/k cq <, 10 4 , but Tcamb /Teh decreases as k increases 
until it reaches a new plateau: Tcamb(^) — 0.77TEH(fc) 
for k/k cq <; 5 x 10 5 . This suppression results from the 
slower growth of these small-scale perturbations after re- 
combination, which is not included in TehCO- 

We base our scale-dependent growth function on the 
ratio Tcamb/Teh at z — 50; for z > 3, modes with 
k/k eq <; 10 5 are proportional to T>(a), while modes 
with k/k eq 10 4 are proportional to (2/3) + (a/a eq ). 
Intermediate scales smoothly interpolate between these 
two values; since we will primarily be concerned with 
modes with k/k cq > 10 7 , the details of the transition 
are unimportant. At later times, the cosmological con- 
stant suppresses perturbation growth, and we assume 
that all scales follow the standard growth function D(a) 
for z ^ 3. An explicit expression for our growth function 
is given in Appendix [Cl 

To modify the transfer function to include a pe- 
riod of scalar domination, we multiply it by the ratio 
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FIG. 7: The matter transfer function in the absence of free- 
streaming effects for four values of the reheat temperature: 
Trh = 8.5 MeV, 85 MeV, 780 MeV, and 6.0 GeV, correspond- 
ing to kjiH /k cq = 10 7 , 10 8 , 10 9 and 10 10 , respectively. Pertur- 
bation modes that enter the horizon while the Universe is 
radiation dominated have T(k) oc ln[fc]/fc 2 , but modes that 
enter the horizon prior to reheating have a constant-valued 
T(k). 



<5dm(TR H )/<5dm(7R H = oo), where S dm is evaluated af- 
ter matter-radiation equality and before baryon decou- 
pling. From Eq. (|20p . we see that, for modes with 
k 3> fc cq , including a period of scalar domination takes 
T{k) -> R{k)T(k), where 



R(k) = 



A(k) In 


f4\W -B(fe)a cq 
\e 3 ) a hor (fc) 


9.11 In 





(27) 



In this expression, ahor(k) is given by Eq. and A(k) 
and B{k) are given by the fitting functions in Appendix 
Of°r k > 0.05fcRH- For k < 0.05fcRH, <5dm is n °t affected 
by the period of scalar domination, as seen in Fig. [5] 
Therefore, we take R(k < 0.05fcRn) = 1- The result- 
ing transfer functions (evaluated at z = 0) are shown in 
Fig. ([7]) for four values of the reheat temperature. As 
expected from Fig. [5J we see that the transfer function 
deviates from its standard value when k exceeds /crh, 
and it is scale- invariant for k IO/crr- 



B. Free-Streaming Effects 

We have thus far assumed that the dark matter has 
no momentum when it is produced from the decay of the 
scalar field. We now relax that assumption and introduce 
(wrh): the average velocity (dx/dr) of the dark matter 
particles when the Universe became radiation dominated. 
We will continue to assume that the dark matter is always 
non-relativistic, so (urh) *C 1. 



If the dark matter particles have momentum, then 
their free-streaming will tend to erase structures on scales 
smaller than the free-streaming horizon Af s h [HE HI, [HJ : 



Afsh(i) 



a 



(28) 



where (v) = (vrh ) ( a RH / a ) is the average velocity of 
the dark matter particles after reheating. The free- 
streaming horizon does not change significantly after 
matter-radiation equality, so we may neglect dark energy 
and take 



Afeh (a) 



da' 



+ a'/a c 



(29) 



where 17 r is the present-day radiation density divided by 
the critical density. This integral can be evaluated ana- 
lytically [63 ]: 



Af s h (a) 



2(^rh)«rh 



sinh 



ORH 



sinh- 1 ,/^ 
a 



(30) 

The effects of free-streaming on the matter power spec- 
trum can be approximated by introducing a Gaussian 
cutoff to the transfer function: 



T{k) = exp 



9k 2 

ZK fsh 



lb (A:), 



(31) 



where Tb(fc) is the transfer function without free- 
streaming calculated in the previous section, and 
fcfsh = Xjj£ [H |H, |H, If kss/km > l, then the 

free-steaming cutoff in the transfer function will exponen- 
tially suppress the scale- invariant portion of T (k) given 
by Eq. (|25|) . The growth of density perturbations during 
scalar domination will be completely erased by the free- 
streaming of the dark matter particles after reheating. 
The ratio fcRH~/fcfsh depends only weakly on the reheat 
temperature; if we neglect changes in the number of rel- 
ativistic species after reheating, 



fcRH 
fcfsh 



2(vrh) 



(32) 

For a wide range of reheating temperatures (10 Mev ^ 
T RH <, 10 GeV), fc RH /fe fsh ~ (w RH )/0.06. 

If the dark matter produced via decay does not inter- 
act with the radiation fluid, then a low (i>rh) will occur 
if the mass of the dark matter (or more generally the to- 
tal mass of the decay products) is nearly degenerate in 
mass with the decaying scalar. If the dark matter instead 
interacts with the plasma through a higher-dimensional 
effective operator (such as a four-fermion operator) , then 
we generically expect a larger interaction cross section 
at higher energies. In this case, we expect dark matter 
particles born with large momentum to rapidly lose en- 
ergy to the plasma until interactions become inefficient 
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and the diffusion length to be far shorter than the free- 
streaming length would be without interactions. Calcu- 
lating the evolution of dark matter perturbations includ- 
ing the thermalization and diffusion of high-energy dark 
matter particles goes beyond the scope of this work, but 
after interactions become inefficient we expect the for- 
malism developed here will describe subsequent pertur- 
bation evolution with (vrh) ~ v -E-kd / "Jdm set by the 
energy at which particles typically kinetically decouple 
from the plasma. Finally, we mention that if some frac- 
tion or all of the dark matter did not originate from scalar 
decay and was present during the scalar-dominated era, 
then these particles likely have low velocity [67]], and we 
expect that they would experience the same perturbation 
growth described at the beginning of Section IIVI 
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V. PRIMORDIAL STRUCTURES 

We have shown that matter perturbations on scales 
smaller than the horizon size at reheating and larger than 
the matter particle's free-streaming horizon are enhanced 
relative to their amplitudes in a universe that was al- 
ways radiation dominated. We will now use the Press- 
Schechter formalism [68[ to analyze the ramifications this 
enhancement has on the formation of dark matter halos. 
We must first compute the rms density perturbation in 
a sphere containing an average mass M: 



a 2 (M,z) = J 



A3h. 

-^{D 5Q (k 7 z)T 50 (k)} 2 P p (k)F 2 (kR), 



(33) 

where Tso(fc) is the transfer function evaluated z = 50; 
D 50 (z, k) is the scale-dependent growth function de- 
scribed in Appendix [C] Pp(k) is the power spec- 
trum of superhorizon density perturbations during 
radiation domination; and F(kR) is a filter func- 
tion that suppresses contributions from modes with 
fc- 1 < R = [3M/(47rp m , )] 1/3 , where p m fl is the present- 
day matter density. We determine P p (k) by setting 
the power spectrum of superhorizon curvature fluc- 
tuations A^(fc) = 2.44 x 1(T 9 x (fc/fco)"- 1 , where 
fco = 0.002 Mpc" 1 and n is the scalar spectral index: 
n = 0.963 ±0.012 [nj. 

It is customary to use a spherical top-hat window func- 
tion in real space as a filter when computing a(M); this 
choice implies F(kR) = 3[sm(kR)-(kR) cos(kR)]/(kR) 3 . 
However, this filter function is problematic when the 
effects of a scalar-dominated era arc considered. We 
showed in Section IIV Al that T(k) is scale invariant for 
k iS fcRH- Since P p (k) oc fc™, the integrand in Eq. (f33|) 
for fc > fcRH is proportional to k n+2 F 2 (kR)Ak. With 
the top-hat filter function, F 2 {kR) ~ cos 2 (kR)/(kR) 4 
for fc ^> so for kR 3> 1 and fc <; fcRH, the integrand 
in Eq. ([33| is proportional to fc"~ 2 dfc. Since n is only 
slightly less than unity, this integral is nearly divergent, 
and the portion of the integral with fc 3> 1/R makes a 
significant contribution to a 2 (M). In short, the standard 



FIG. 8: The present-day rms density perturbation a in a 
sphere containing an average mass M for three values of the 
reheat temperature: Trh = 8.5 MeV, 85 MeV and 780 MeV. 
For these values of T RH , Mrh = 54 M e , 0.054 M e , and 
5.4 x 10 _5 M®. We also show cr(M) for reheat tempera- 
tures greater than 12 GeV; in the mass range shown here, 
the presence of a scalar-dominated era does not affect a(M) 
if Trh > 12 GeV. Note that a(M <; M RH ) depends only 
weakly on M and is independent of Trh, while a(M ^ Mrh) 
is significantly enhanced compared to a for higher values of 
T RH . Also, for all values of Trh, a{M & Mrh) oc M- (n+3)/6 , 
where n = 0.963 is the scalar spectral index. 



top-hat filter fails to prevent modes with fc -1 <C R from 
contributing to a 2 (M) when T(k) is scale- invariant for 
large fc. 

To solve this problem, we use a different filter function 
when calculating a 2 (M): we convolve the spherical top- 
hat window function with a Gaussian window function 
with a much smaller radius. The resulting filter function 
is 



F(kR) 



exp 



(kR) 3 



-k 2 (aRf 



sm(kR) - (kR) cos(fci?)] (34) 



with a <C 1. For a ^ 10 -4 , this window function is 
virtually indistinguishable from a top-hat of radius R in 
real space, and a(M) for R » 1/A;rh changes by less than 
0.01% compared to its value with a = 0. This robustness 
prompts us to set a — 10 -4 when computing a(M). We 
note, however, that the variance on smaller scales is more 
sensitive to the choice for a; for R -C 1/fcRH, changing 
a by a factor of ten in either direction changes a(M) by 
about ten percent. 

The filter function ensures that a(M) is sensitive to 
scales that enter the horizon during scalar domination 
only if R(M) <, k^. Consequently, it is useful to define 
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FIG. 9: The rms density perturbation a in a sphere containing 
an average mass M for four values of the mean dark matter 
particle velocity at reheating (urh), in units where c = 1. The 
reheat temperature is 8.5 MeV. If the dark matter particles 
have a nonzero velocity dispersion, free-streaming erases den- 
sity perturbations on scales smaller than the free-streaming 
horizon, making a(M) nearly constant for small masses. 



a characteristic mass Mrh such that R(Mnn) = fc R ^ 
1 0MeV \ 3 / 9 *sPrh1\ / 10.75 

t rh J V 



M RH = 32.7A/ a 



10.75 



g* Prh] 



3/2 



(35) 

For M > AfRH, the filter function effectively restricts the 
integral in Eq. (f3"3"|) to k <5 fcRH, and cr(M) is insensitive 
to the reheat temperature. In contrast, if M < Mrh, 
<t(M) is most sensitive to scales with k > £rh- The 
transfer function is scale-invariant at these scales, so the 
density power spectrum is a simple power law P(k) on k n . 
It follows that cr(M £ Mrh) oc M~ (n+3)/6 , as shown in 
Fig. E 

As discussed in Section IIV B| we can include the ef- 
fects of the dark matter particles' random motions by 
adding a Gaussian cutoff to the transfer function, as in 
Eq. pip. Since the resulting transfer function is no longer 
scale-invariant for k > fcf s h, we expect a(M <S Mrh) to 
deviate from the power law behavior shown in Fig. [8] 
when M <, M fsh , where R(M ish ) = k~^ = X ish . Fig- 
ure [3] shows that this is the case; in this figure, Trh is 
8.5 MeV, and er(Af) is shown for four values of (urh)- 
For (u RH ) = 0.01,0.001, and 0.0001, Eq. ([30]) implies 
M fsh = 0.31 M m , 3.1 x 10" 4 M ffi , and 3.1 x 10" 7 M e , re- 
spectively, and we see that a(M) is nearly constant for 
M <> Mf s h. Figure [S] also shows that the growth of the 
matter density perturbations prior to reheating is com- 
pletely erased by free-streaming if (wrh) ^ 0.01; this is 
also true for other values of Trh- 

We have seen that the growth of matter perturbations 
during the scalar-dominated era leads to an enhancement 
cr(M) for M <, Mrh, provided that (v RH ) <, 0.01. Since 
Mrh ^ 1000 M®, the amplitude of a(M < Mrh) de- 
pends on the amplitude of the scalar power spectrum 



on scales that are much smaller than the pivot scale 
fc = 0.002 Mpc" 1 . Consequently, a(M <, Mrh) is very 
sensitive to changes in the scalar spectral index n. We 
can quantify the dependence of cr(M) on n by noting that 



g(M,m) 
<j(M, n 2 ) 



1 



k R a {M) 



(36) 



where R a (M) is the average scale that contributes sig- 
nificantly to cr(M). For M ^ M RH , R a {M) = R(M) 
because the transfer function and filter function sup- 
press the contribution from smaller scales. For M <, 
Mrh and (i>rh) = 0, all scales between aR and R 
contribute equally, so R a {M) ~ R(M)y/a. Free- 
streaming introduces an additional small-scale cutoff; for 
aR & A fsh <, R, R a (M) ~ y/R{M)X bh , and R a (M) ~ 
Af s h if R ^ Af s h- For example, if Trh = 8.5 MeV and 
(vrh) = 0, ct(0.001 Mrh) increases by 37% if n is 0.987 
instead of 0.963. If we instead take (wrh) = 0.001, then 
<r(0.001 A/rh) increases by 32% for the same change in 
n; the introduction of a free-streaming cutoff has reduced 
effects of increasing the scalar spectral index. 

Now that we have calculated cr(Af), we can use the 
Press-Schechter formalism 68] to obtain a halo mass 
function: 



dn 



din M 



2 Pm,0 
7T ~W 



din a 



d In M 



a(M,z 



■ exp 



2cr 2 (M, z) 



(37) 

where n is the comoving number density of halos with 
mass M at redshift z, and S c is the critical linear over- 
density. For z ^ 2, S c = 1.686, and it decreases slightly 
at smaller redshift. Since an era of scalar domination 
changes a(M) for M <, Mrh, only the abundance of 
these small-mass halos will be affected. These small-mass 
halos form at very high redshift, and then they merge to 
form larger-mass halos. The Press-Schechter mass func- 
tion does not account for subhalos, however, so it will be 
most useful for us to examine the halo mass function at 
high redshift, before halos with M <, Mrh are absorbed 
into larger halos. Rather than consider the number den- 
sity of these objects, it is more enlightening to compute 
the fraction of the mass that is contained in these objects 
as a function of redshift: 



df 



M dn 



d\nM p m , d\nM 



(38) 



The ratio S c /<j plays a pivotal role in the Press- 
Schechter mass function, and it is useful to define a func- 
tion M*(z) such that cr(M*,z) = 5 C . Figure [TU] shows 
M*(z) for different values of the reheat temperature and 
(urh)- We see that M*(z) decreases precipitously with 
increasing redshift for M <; Mrh. For smaller masses, 
M*(z) transitions to a power law. Free-streaming by the 
dark matter particles steepens this power law slightly be- 
cause a(M) increases less quickly with decreasing M if 
(vrh) is nonzero, as shown in Fig. [5J Also recall that 
cr(M) is nearly constant for M <5 Mf s h; this constancy 
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FIG. 10: The value of M* such that cr(M«) = 5 C plotted as a 
function of redshift for different values of the reheat tempera- 
ture Trh and the average particle velocity at reheating, (wrh), 
in units where c = 1. For T RH = 8.5 MeV and T RH = 85 MeV, 
the reheat horizon mass is Mrh = 54 M® and 0.054 M®, re- 
spectively. We see that M*(z) is a power law for M* & Mrh- 



would cause M*(z) to decrease rapidly with increasing 
redshift if M* ^ Mf s h- We see in Fig. [TU1 however, that 
if («rh) £ 0.001, M* £ M fsh for all z < 500. There- 
fore, we expect that free-streaming will not prevent struc- 
tures from growing at redshifts z <, 500, provided that 

(«BH>£ 0-001. 

Before we consider the effects of free-streaming on the 
halo mass function further, we examine the mass func- 
tion with (urh) = 0. If a(M) is a power law, as it 
is for M <C Mrh, then the mass function depends on 
redshift only through the ratio M/M*(z). In this case, 
the differential bound mass fraction df/dhxM peaks at 
M = M*; for dlna/dlnM = — (n + 3)/6, the maximum 
value is df /d\nM\M, = 0.32, and 31% of the dark mat- 
ter is contained in halos with M* < M < 5M*. Since 
the abundance of halos with M <; M* is exponentially 
suppressed, a negligible fraction of the mass is contained 
in halos with M > 5M*. Provided that 5M* < M RH , 
the fact that a(M) is not a power law for M <; Mrh will 
not affect the bound fraction. 

The 5Af* -C Mrh case is exemplified by the z = 100 
curve in Fig. 111! which shows the differential bound 
fraction df/dlnM for T RH = 8.5 MeV. For this reheat 
temperature, Mrh = 54M®, and Fig. [TU] shows that 
M* = 0.18M e at z = 100. Integrating the z = 100 curve 
in Fig. [TT] reveals that almost half of the dark matter is 
contained in halos with M > O.IM® at this redshift. In 
the standard cosmology, the fraction of the dark matter 
contained in such halos is only 10~ 10 at z = 100. As the 
redshift increases past 100, the differential bound frac- 
tion will keep the same shape as the z = 100 curve while 
the peak slides to smaller masses. The other curves in 
Fig. [TT] show what happens to the bound mass function 
as the redshift decreases. While z ^ 11, M* > Mrhi 



FIG. 11: The differential fraction of the dark matter mass 
that is bound into halos of mass M [see Eq. (|38f) ] plotted 
at several redshifts for Trh = 8.5 MeV. While z <; 11, the 
critical mass M» is smaller than the reheating horizon mass 
Mrh, and the peak moves to larger masses as the redshift 
decreases. When z < 10, M* > Mrh, and the peak remains 
fixed at 2M®. The peak decreases in amplitude because the 
fraction of mass contained in these microhalos decreases as 
they are absorbed into larger halos. 
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FIG. 12: The differential fraction of the dark matter mass 
that is bound into halos of mass M plotted at several redshifts. 
The solid curves show the bound fraction for Trh = 8.5 MeV, 
while the dashed curves show the bound fraction in the ab- 
sence of a scalar-dominated era. The two scenarios are indis- 
tinguishable for M <; 1000 M® . The abundance of halos is 
higher for halos with M 1O 1O M0 because baryons fall into 
these halos, making them grow faster. 



and the peak in df/dlnM follows M*, moving to larger 
masses as the redshift decreases. The peak height de- 
creases because \d\na/d\n M\m, is decreasing, as seen 
in Fig. HI 

When z £ 11, M* > M RH for T RH = 8.5 MeV, 
and larger-mass halos begin to form. These larger ha- 
los absorb some of the microhalos formed at higher red- 
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shift, and so the fraction of mass in these microhalos 
decreases, as seen in Fig. [TT] The transfer of mass from 
small halos to larger halos is more apparent in Fig. IT2| 
which also shows the differential bound fraction without 
a scalar-dominated era. We see that a period of scalar 
domination does not affect the abundance of halos with 
M ^> Mrh- Without a scalar-dominated era, however, 
halos with M Mrh are far less abundant and they 
form much later. For example, if Trh = 8.5 MeV, the 
abundance of Earth-mass microhalos peaks at a redshift 
of z = 32 when df /d\nM\ MfSj — 0.25. In contrast, iso- 
lated Earth-mass microhalos are most abundant at a red- 
shift of z = 12 in the standard scenario, and even then 
df/d\nM\u @ = 0.010. 

Given our understanding of the bound mass fraction 
when (urh) = 0, we can now examine the effects of 
free-streaming. Figure [TU] shows that free-streaming de- 
creases M* when M* <, Mrh, and Fig. [S] shows that 
free-streaming decreases d In a/d In M when M £ Mrh- 
We conclude that free-streaming will shift the peak in the 
bound mass fraction to smaller halo masses and decrease 
its amplitude. Figure [T3l confirms that this is indeed the 
case. We see that free-streaming delays the formation of 
microhalos of a given mass and suppresses their abun- 
dance; for example, if (vrh) = 0.001 and Trh = 8.5 
MeV, the abundance of Earth-mass microhalos peaks at 
a redshift of z = 22 when df/d\nM\ MfB = 0.17. This is 
still significantly higher than the standard abundance of 
these microhalos. We conclude that free-streaming does 
not erase the structural imprints of an extended scalar- 
dominated era, provided that (wrh) ^ 0.001. Microhalos 
with M < A^rh still form earlier and in far greater num- 
bers than they would in the standard scenario, but free- 
streaming does slightly delay their formation and sup- 
press their abundances. 



VI. SUMMARY AND DISCUSSION 

The content of the Universe prior to the onset of big 
bang nucleosynthesis is unknown. It is possible that the 
inflaton continued to dominate the energy density of the 
Universe long after inflation ended. After inflation, the 
inflaton oscillates around the minimum of its potential, 
and if the potential is quadratic around its minimum, 
then the time-averaged pressure of the scalar field van- 
ishes. It is also possible that the inflaton decayed quickly 
to radiation, and then a second oscillating scalar field 
may have come to dominate the Universe, as in some 
versions of the curvaton scenario [23j-|26j and in several 
realizations of string theory |20l-t22|. Finally, the infla- 
ton may have decayed into a short-lived heavy particle 
that dominated the energy density of the Universe prior 
to its decay. In all of these scenarios, the Universe was 
effectively matter dominated prior to nucleosynthesis. 

During this early matter-dominated era, subhorizon 
perturbations in the dominant component of the Universe 
grow linearly with the scale factor [39M42| . We have in- 
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FIG. 13: The differential fraction of the dark matter mass that 
is bound into halos of mass M for different values of the mean 
dark matter particle velocity at reheating, as a fraction of the 
speed of light (urh). The reheat temperature is Trh = 8.5 
MeV, and the function is shown at two redshifts. We see that 
free-streaming moves the peak in the bound mass fraction to 
smaller masses and decreases its amplitude. 

vestigated what happens to these density perturbations 
after the dominant component during the early matter- 
dominated era decays and the Universe becomes radia- 
tion dominated. We model the reheating of the Universe 
using three fluid components: a pressureless "scalar" that 
decays into a tightly-coupled radiation plasma and de- 
coupled nonrelativistic dark matter. These dark matter 
particles may have a small, but nonzero, velocity disper- 
sion when they are created, and we examine the effects 
of free-streaming on their perturbations. 

We find that, although subhorizon radiation perturba- 
tions grow during the scalar-dominated era due to the 
coupling between the radiation and the scalar field, this 
growth is erased after the Universe becomes radiation 
dominated. When the scalar field energy density van- 
ishes, these radiation perturbations begin to oscillate, 
but the amplitude of these oscillations is much smaller 
than the maximum value of the radiation density per- 
turbation during the scalar-dominated era. Moreover, 
the amplitude of the oscillations in the radiation density 
perturbations for modes that enter the horizon well be- 
fore reheating is significantly smaller than the amplitude 
of modes that enter the horizon after the Universe be- 
comes radiation dominated. An early era of matter dom- 
ination suppresses perturbations in the radiation density 
field on scales smaller than the horizon scale at reheat- 
ing. If the dominant form of dark matter is not cre- 
ated directly from reheating but instead from the radia- 
tion bath after reheating, this radiation suppression re- 
sults in a suppression of the matter power spectrum for 
k > 23 A (Trh/2 MeV) kpc -1 . Depending on the details 
of the dark-matter microphysics, this reheating suppres- 
sion may be the primary cutoff in the matter power spec- 
trum. This could include corners of WIMP parameter 
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space with inefficient kinetic coupling [56| , FIMP models 
[57j . and hylogenic models (121 ] . 

The effects of reheating on dark matter produced di- 
rectly from the decay of the scalar can be much more 
dramatic. Unlike the radiation density perturbation, the 
perturbation in this dark matter density is enhanced by 
the growth it experiences during the scalar-dominated 
era. While the Universe is scalar dominated, dark matter- 
density perturbations grow linearly with the scale factor 
after they enter the Hubble horizon. When the Universe 
becomes radiation dominated, the dark matter pertur- 
bations grow logarithmically from their amplitude at the 
moment of reheating. Consequently, perturbations on 
scales that enter the horizon prior to reheating are sig- 
nificantly enhanced compared to larger-scale perturba- 
tions, and the matter power spectrum on small scales is 
proportional to k n , where n is the scalar spectral index. 
Free-streaming breaks this scale invariance and exponen- 
tially suppresses the power spectrum on scales smaller 
than the free-streaming horizon. We found that free- 
streaming completely erases the perturbations that grow 
prior to reheating if the average velocity dispersion of the 
dark matter particles at reheating is greater than 0.01c. 

We used the Press-Schechter halo mass function [68[ 
to investigate how the growth of perturbations prior to 
reheating affects the formation of dark matter halos. As 
expected, only halos with masses less than the horizon 
mass at reheating are impacted. The reheating horizon 
mass (Mrh) is proportional to Xrij and equals 260Af ffi 
if the reheat temperature was Trh = 5 MeV. We found 
that microhalos with M < Mrh form much earlier than 
they would if the reheat temperature were higher, and 
they contain a much larger fraction of the dark mat- 
ter. For instance, if Trh = 8.5 MeV, then the Press- 
Schechter mass function predicts that 15%, 43% and 65% 
of the dark matter is contained in microhalos with masses 
greater than 1O _6 M0 at redshifts of 100, 50, and 25, re- 
spectively. In contrast, the corresponding fractions are 
10~ 10 , 0.04% and 5% if the reheat temperature is higher 
than 100 MeV. Numerical simulations of the formation 
of the first microhalos confirm this prediction; they find 
that only 1.5% of the dark matter is in bound in these 
microhalos at a redshift of 31 (6{|, and 5% of the dark 
matter is in these microhalos at a redshift of 26 [7C§ . 

We conclude that low-temperature reheating with di- 
rect dark matter production results in an abundance of 
earth-mass or smaller microhalos. Prior to the formation 
of larger dark matter halos, these microhalos contain a 
significant fraction of the dark matter. What happens to 
these microhalos as they merge to form larger halos, and 
do they survive until the present day in the Milky Way's 
halo? The fate of the earth-mass microhalos that form in 
the standard cosmological scenario has been studied ex- 
tensively d^-lzl] ; these subhalos survive their absorption 
into larger halos, but then they lose most of their mass 
due to interactions with stars. The survival probability 
for a microhalo near the Sun has been calculated to be as 
high as 0.17 [72j], but this factor depends strongly on the 



microhalo's orbit, and some studies predict that nearly 
all the local microhalos are either destroyed by stellar en- 
counters (Til [73|, [75| or lose a significant fraction of their 
mass [74], |76| . In the low-temperature reheating scenarios 
considered here, however, these microhalos are far more 
numerous than in the standard scenario; even if there 
is a high probability that they do not survive inside the 
Galaxy, there may still be a sizable population in our 
stellar neighborhood. 

Furthermore, the microhalos in the low-temperature- 
reheating scenario form far earlier than their standard 
counterparts, and consequently, they are far denser. The 
virial density of a microhalo is proportional to the critical 
density at the time of its formation; at high redshifts, the 
mean density of a microhalo that formed at a redshift z f 
is therefore proportional to (1 + z/) 3 - Ref. found 
that the microhalo survival probability near the Sun in- 
creases sharply as the mean density within the clump 
increases; while a microhalo that formed at a redshift of 
65 has a survival probability of 0.17 at the Sun's loca- 
tion, a microhalo that formed at a redshift of 120 has 
survival probability of 0.5, and a microhalo that formed 
at a redshift of 300 has survival probability of 0.9. It is 
therefore reasonable to expect that many of the microha- 
los produced in the low-temperature-reheating scenario 
are intact subhalos today. 

How could we detect these small and dense subha- 
los? While a complete analysis of their observational 
signatures lies beyond the scope of this work, we briefly 
highlight a few promising detection avenues here. If 
the dark matter self-annihilates, these compact micro- 
halos will be gamma-ray sources; they may be ob- 
served as point sources, and they would make a sig- 
nificant contribution to the observed gamma-ray back- 
ground [H, [zO, [Z3, [ll, ■ Unfortunately, since the 
dark matter particles in these microhalos were created 
through scalar decay and cannot be thermal relics, there 
is no lower limit on their self-annihilation cross section. 
Numerous and dense subhalos near our location would 
also affect direct detection rates by changing the local 
dark matter density [76|, [8l[ and the direction of the dark 
matter wind which would alter both the event rate 
and its annual modulation. If dark matter particles are 
directly detected, then the temporal correlation of the 
signal can be used to probe the abundance of local small 
subhalos because the event rate will change as a micro- 
halo passes through the Solar System [80j. 

Of course, it is entirely possible that the dark matter 
particles in these microhalos do not self- annihilate and in- 
teract too weakly to be detected directly. The only guar- 
anteed observational signatures of dark matter microha- 
los are gravitational. Dark matter microhalos are grav- 
itational lenses; we could hope to detect them through 
their impact on individual images of a strongly lensed 
quasar [7(1 Isl, [83| . Unfortunately, even though the early- 
forming microhalos that result from perturbation growth 
during reheating are more compact than their standard 
counterparts, their virial radii are still much larger than 
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their Einstein radii, and it is unlikely that they could 
act as strong lenses. Recently, it has suggested that we 
may be able to detect subhalos within our galaxy through 
their astrometric microlensing signatures [84j or their im- 
pact on the pulse arrival times from millisecond pulsars 
[85j- These studies found that subhalos in the standard 
scenario are too rare and too diffuse to be detectable 
through these methods given our current level of astro- 
metric and timing precision, but the abundant and dense 
microhalo population resulting from low-temperature re- 
heating with direct dark matter production would signifi- 
cantly enhance both signals. It is also possible that these 
microhalos may be dense enough to be detected through 
photometric microlensing (86[. 

Finally, we note that the Laser Interferometer Space 
Antenna (LISA) will be able to detect compact objects 
with masses greater than 10 16 g that pass near its de- 
tectors; if these objects are the primary component of 
the dark matter, then a few such events are expected per 
decade [87|- LISA may also be sensitive to small and 
dense dark matter microhalos, if their radii are smaller 
than LISA'S arm length (~5x 10 11 cm) %7%. Microhalos 
that form before z = 340 can have virial masses greater 
than 10 16 g and virial radii less than 5 x 10 11 cm. Such 
microhalos are less massive than the earth-mass micro- 
halos that we have focused on and would form at high 
redshift for larger values of the reheat temperature. For 
instance, if Trh = 12 GeV, then the critical mass at 
z = 500 is M* = 2.5 x 10 16 g, which is much less than 
the reheat horizon mass (Mrh = 4 x 10 19 g). The Press- 
Schechter mass function predicts that 30% of the dark 
matter is contained in halos with M ~ M* if M* <C -Mrh- 
If these microhalos survive their absorption into larger 
halos, they could make up a significant fraction of the 
Galaxy's dark matter halo, and they could be detectable 
by LISA. It would be interesting to explore this possibil- 
ity, and the possibility of detecting early-forming micro- 
halos through astrometric, photometric and timing mi- 
crolensing in more detail; detecting these small structures 
or constraining their abundance would probe the origins 
of dark matter and the state of the Universe between 
inflation and nucleosynthesis. 
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Appendix A: Derivation of perturbation equations 

We treat the oscillating scalar field, the radiation, and 
the dark matter as perfect fluids with energy momentum 
tensors 



= (p + p)u»u v +pff 



)J,V 



(Al) 



where p and p are the fluid's density and pressure, re- 
spectively, and = dx^/dX is its four- velocity. In this 
Appendix, we will take x° to be the proper time t. The 
dark matter and the oscillating scalar fields are both pres- 
sureless fluids, while the radiation has p — p/3. Since 
the scalar field is decaying into radiation and matter, 
these three fluids are exchanging energy, as described in 
Eqs.fl])-®. We can express this energy exchange covari- 
antly: 



\ ( (i) T^) =Q« 



where i denotes the individual fluids. 
Eqs.CQ)-© that 



(A2) 

It follows from 

(A3a) 
(A3b) 
(A3c) 



Thus we see that, in our three-fluid model, 



(dm) 



0, 



(A4) 



as required by the conservation of energy and momen- 
tum. 

We obtain the perturbation equations by evaluating 
Eq. (|A2p with the perturbed metric 

ds 2 = -(1 + 2*)di 2 + a 2 (i)<%(l + 2$)dx i dx j (A5) 

and with perturbations in the density of each fluid: 
Pi(t,x) = pj(t)[l + Si(t,x)]. We also introduce pertur- 
bations to the four-velocity of each fluid: u° = (1 — \P) 



and 
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(1 



where VjL. = dx J /dt is the fluid 



velocity of the ith fluid. It follows that 



(A6) 
(A7) 



to first order in the perturbations. Thus we see that 
is a first-order quantity, while has both a zero-order 
r^pj] and a first-order component 



W),(o) _ 



component [Qq 

[Q^ ),{1) =i>^ + *)]. 

The p = component of Eq. 
each fluid 

dS , J , N d$ 
— + (l + w)- + 3(l + w) — 
dt a dt 



2|) implies that, for 
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where w = p/p is the fluid's equation of state parame- 
ter, 9 = adiV % is the divergence of the fluid's conformal 
velocity, and Q^ 1 and Q^P are the zero-order and first- 
order components of Qa for this fluid. The divergence of 
the spatial components of Eq. (IA2|) implies that 



V 2 <5 



^- + (l-'3w)H8+ — + - 

dr a 1 + w a 



1 



o(l + iu) 



Qo 0) * 



(A9) 



Evaluating these equations for each fluid using Eq. (|A3 
yields the perturbation equations given by Eq. (|10|) . 



Appendix B: Initial conditions for the Perturbations 

The equation suite that we solve is 



a 2 E{a)8' lt> {a) + 9^{a) + 3a 2 E(a)&{a) 
a 2 E (0)6^(0) + aE(a)9^ + fc 2 $(a) 

a 2 E(a)S' r {a) + \e r {a) + 4a 2 E(a)& (a) 

a 2 E(a)9' r (a) + ~k 2 ^(a)-k 2 ^- 

a 2 E(a)S' dm (a) + 9 dm (a) + 3a 2 E(a)& (a) 

0^(0)^(0) + a£(a)0 dm + fc 2 $(a) 
fc 2 $ + 3aE 2 (a) [a 2 $» + a$(a)] 



ar $(a), 
0, 



#(«) f 



Pr\ a ) 



3 - 

-0<j,{a) - 8 r (a) 



p°(a) 

/ ~o , \ «r [(50(a) - <5 dm (a) - $(a)] , 



/ -0 7 ^ aT 4 



9</,(a) - 9 dm (a) 



-a 2 [p°(a)54a) + p° r (a)5 r {a) + ^(o)] , 



(Bla) 
(Bib) 

(Blc) 
(Bid) 
(Ble) 
(Blf) 
(Big) 



r 



where a prime denotes differentiation with respect to a. 

We chose the initial time to such that all modes of 
interest are super- Hubble (k < 1). In conformal New- 
tonian gauge, super-Hubble perturbations do not evolve; 
to zeroth order in k 2 we may set = 8'^ = 0. By 

choosing <C 1, we have ensured that the Universe 
is initially dominated by the scalar field energy density, 
which implies that E(a) = aT^I 2 and p$ 3> p r , Pdm- Also, 



P<t>{h 



1 implies that , 



prior to the decay of the 



scalar field. With these conditions, and neglecting terms 
that are proportional to fc 2 , and Eq. ( |Blg[ ) implies 
that S^fto) = 2$ (to)- Equation (jBlbj) reduces to 



v o, 



/c 2 $ 



0. 



(B2) 



The solution to this equation is 6$ = — (2/3)fc 2 $ v / a. 

The solutions for pdm(t) and p r (t) prior to the decay 
of the scalar field [Eqs. (|6]) and (0] imply 



(1 



PrW 

/ -0 ? -, aT 4 
Pdm ( a ) 



2ja' 
3 



2y/E' 



(B3) 



(B4) 



Therefore, the right-hand sides of Eqs. (|Blc|) and (|Blc|l 
are not negligible at t = to, even though <C 1. The 



left-hand sides of Eqs. (jBlcl) and (|Blc|) are initially zero, 
however, because the density perturbations and $ do not 
evolve while the mode is superhorizon, and the velocity 
terms are proportional to k 2 . To solve these equations 
at the initial time, we must set the right-hand sides of 
Eqs. (jBlcp and (jBlel) to zero by demanding that 5 r — 
fidm = ^0 — Therefore, our first-order initial conditions 
for the matter and radiation perturbations are S r (to) = 

^dm (to) = $(*o). 

For superhorizon modes in a universe dominated by an 
oscillating scalar field, Eq. (IBldl) becomes 



2y^ 



= 0; 



(B5) 



we have used 9^ = — (2/3)k 2 $y/a and S r — <f> to ar- 
rive at this equation. The solution to this equation is 
9 r = — (2/3)fc 2< i>v / a. For superhorizon modes in a uni- 
verse dominated by an oscillating scalar field, Eq. (|Blfj) 
also becomes 



Va6' dm (a) +2fc 2 $ 



2v^ 



6a 



0, 



(B6) 



where we have again used 6$ = — (2/3)k 2 $>y/a. Since this 
is the same equation as we obtained for 9 ri we have the 
same initial condition for 9 dm . 
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We can obtain more accurate initial conditions for the 
density perturbations by inserting the early-time solu- 
tions for the velocity perturbations into Eq. (|B1[) . To- 
gether, Eqs. fBlgP and (TBI aft imply that the 0(k 2 ) term 
in the early-time solution for $(a) is proportional to 
a~ 5 / 2 and has an undetermined coefficient. Since this 
is a decaying mode that diverges as a approaches zero, 
we choose this coefficient to be zero and continue to set 
$ = $o and &'(a) = at early times. If we do not neglect 
the fc 2 $ term in Eq. QBlgD , we see that 



2 ~ 

64,(0) = 2$o + ~fc 2 $oa 



(B7) 



at early times. Using this solution and 

r = Odm = -(2/3)P$Va in Eqs. pTc) and (jBTc|) 
gives 



46 - 

6 r (a) = $ + — fc 2 $ a 
63 

2 ~ 

6dm(a) = $0 + gfc 2 $oa 



(B8) 
(B9) 



at early times. Like our initial condition for <£>, these solu- 
tions also include undetermined decaying modes, and we 
set these modes to zero. As long as the Universe is scalar- 
dominated, which implies T^a <C 1, these early-time so- 
lutions for the perturbations satisfy Eq. (IB1[) . with one 
exception: there is a neglected 0(k 4 a) term in Eq. (|Bld|) . 
Therefore, we expect that the radiation perturbation will 
deviate from this early-time solution when the mode en- 
ters the horizon. The scalar and matter perturbations, 
however, will follow these solutions until the Universe 
becomes radiation dominated. 



Appendix C: Fitting Functions for the Transfer 
Function 

We want to find functions A(k/k^) and B(k/k R n) 
such that 



<5 dm (a, k) = y®° 



A In 



Ba 

flhor 



(CI) 



where k = ah OI H(ai lor ). We know that A = 9.11 and 
B = 0.594 for fc/fc RH ^ 1. We also know that A and B 
are given by Eq. (fl9| for /c/Zcrh *5 10. We found that the 
numerical solution for <5dm(a, k) for 0.05 ^5 fc/fcRH ^ 100 
is well-fit by the following expressions for A(k/kR,n) and 
B(k/k RH ): 



A(x) 



exp 



0.609 



{l + 2.15(lna;- 1.52) 2 } 



21 1.38 



(C2) 



9.11 5(5.02 - x) + -x 2 S(x - 5.02) 
5 



lnB(x) = ln(0.594) 5(5.02 - x) + In 



S(x - 5.02), 



where 



S(y) = - 2 



tanh ( — 



1 



(C3) 



serves as a smooth step function. We have tested these 
fitting functions for several values of £;rh and a, and we 
find that Eq. (|C1I) fits the numerical solution for <5dm(a, k) 
to within 5% for 0.05 < k/k^, provided that the modes 
have entered the Hubble horizon and 5p r ^> <5pdm- 

We also want to find a scale-dependent growth function 
D(k, z) such that 



6d m (k, z) = D 5Q (k, z) x 6 dm (k, z = 50). 



(C4) 



The function D$o(k,z) should mimic the scale depen- 
dence of Xcamb/Ieh a ^ re dshift z = 50. This ratio is 
approximately unity for k/k cq <5 10 4 , and then it de- 
creases to 0.773 for k/k cq <; 10 5 . This transition follows 
Tcamb/Teh ^ D s (k), where 



D s (k) = 



D a -D b 



1 



48500 



2.1 



D h 



(C5) 



with D a = 1 and D b = 0.773. For redshifts 
3 ^ z <, 500, the modes with Tcamb/?eh — 1 are pro- 
portional to (2/3) + (a/a oq ), while the modes with 
2camb/Teh — 0.773 are proportional to V(a) as de- 
fined by Eq. (|2Tj) . Therefore, we define our scale- 
dependent growth function in this redshift range as 
D 5Q (k,z) — D s (k) with 



D a (z) 
D b (z) 



2 1 1 + Zpq 

3 ' 1 + z 

2 1 1 + Zcq 

3 ' 51 



V(z) 



V(z = 50)' 



(C6) 
(C7) 



At lower redshifts, the dark energy slows the growth of 
density perturbations. Since we expect the microhalos 
described in Section [V] to be contained within larger 
structures at these redshifts, the transfer function at 
z <5 3 is not very important to our analysis. To facilitate 
comparisons with the standard cosmological scenarios, 
however, we do report T(k) and a(M) evaluated at the 
present day, and so we need to extend our growth func- 
tion to z = 0. We take 



D 50 (k,z) 



D s (k,z) for z > 2.67 

Dl{z) D s (k, z = 2.67) for z < 2.67 

(C8) 



£>i(z=2.67) 

where D s (k z ) is defined by Eqs. JC5])-([C7]), and 
5n M H(z) rV+'r 1 



Di(z) 



H 3 

H ° da. 



2 H J [aH(a)}3 
We chose z = 2.67 to be the transition point because 
D t (z = 2.67) 



(C9) 



D x (z = 50) 



D a (^ = 2.67). 



(C10) 
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Consequently, on large scales (k/k eq & 10 4 ), 

^o(,<2.67) = ^L y . (Oil) 
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